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ABSTRACT 

This study investigates the expansion properties of nickel (Ni) bubbles in 
Type la supernovae (SNe la) due to radioactive heating from the ^^Ni 
^^Co — i> ^^Fe decay sequence, under the spherically symmetric approximation. 
An exponentially-declining medium is considered as the ejecta substrate, allowing 
for the density gradient expected in a Type la supernova (SN la). The heating 
gives rise to an inflated Ni bubble, which induces a forward shock that com- 
presses the outer ambient gas into a shell. As the heating saturates, the flow 
tends toward a freely-expanding state with the structure frozen into the ejecta. 
The thickness of the shell is ~ 1% of the radius of the bubble, and the density 
contrast across the shell reaches x ^ 100 in a narrow region that is limited by 
numerical resolution. 

In the adiabatic case (in which no radiative energy diffuses across the shell) the 
structure of the shell can be approximately described by a self-similar solution 
that is determined by the expansion rate and ambient density gradient of the 
shell. In the radiative case, the shell expansion weakens but remains comparable 
to the adiabatic solution. The density contrast of the inferred ejecta clumps (shell 
components) is greater than that given by the model that uses a uniform ejecta 
substrate, while the interaction of the clumps with the remnant is delayed to a 
more advanced evolutionary stage. The explosion parameters of the SN are varied 
to examine whether the created clump characteristics are consistent with those 
of the ejecta knots that are present at the edges of Tycho's remnant. Explosion 
conditions similar to those of successful explosion models are found, including 
the deflagration W7 and the delayed detonation yield with sufficient ejection 
velocities as well as timely initiations for the clump-remnant interaction, whereas 
the luminous helium (He) detonations and the low-energy Chandrasekhar-mass 
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explosions are unfavorable. The clump speed can be increased as the initial 
density contrast of ^^Ni, 0.5 < a; < 1, is reduced, as determined by a realistic 
elemental distribution. In all cases, the occurrence of the reverse shock (RS) 
impact of the clump with the remnant is expected in under 2000 yrs after the 
supernova explodes. 

Subject headings: hydrodynamics - nuclear reactions, nucleosynthesis, abun- 
dances — supernova remnants — supernovae: general 

1. INTRODUCTION 

Observations of SN 1987A have shown that the distribution of Fe in the ejecta is not 
what would be expected in the simplest models; it has higher velocities than expected and 
a large filling factor for its mass of 0.07 M©, as determined from the supernova light curve 
(McCray 1993; Li et al. 1993). A plausible mechanism for the large filling factor is the 
Ni bubble effect, in which the radioactive progenitors of the Fe expand relative to their 
surroundings because of the deposition of radioactive power (Wooslcy 1988; Li et al. 1993; 
Basko 1994). This effect is important during the first ~ ten days after the supernova is 
formed, when the radioactive power is significant and the diffusion of energy has not yet 
become important. 

An expected effect of the Ni bubble expansion is to create clumps in the nonradioactive 
gas of intermediate elements. Widespread evidence indicates the cjccta of core-collapse 
supernovae are clumpy. The oxygen line profiles in the nearby Type 11 supernovae SN 
1987A (Stathakis et al. 1991) and SN 1993J (Spyromilio 1994; Matheson et al. 2000) 
showed evidence of the structure, implying that the gas is clumped. The velocity range of 
the emission extends to 1, 500 km s"^ in SN 1987A and 4, 000 km s'^ in SN 1993 J. Clumping 
was also evident in the prototypes of the oxygen-rich supernova remnants (SNRs) Cassiopeia 
A (Cas A) (Hughes et al. 2000; Hwang et al. 2000) and Puppis A (Winkler et al. 1988), 
and in the Type lb SN 1985F (Filippenko & Sargent 1989). As for SNe la, observations of 
SN1006 suggested that Si freely expands in the velocity range 5,600 - 7,000 km s'^ while 
its Fe mass is only estimated as 0.075 — 0.16 Mq (Hamilton et al. 1997). Moreover, in 
Tycho's SNR (SN 1572), X-ray observations revealed two knots, one enriched in Si and the 
other in Fe, protruding non-deceleratingly from the SE edge of the remnant at a velocity of 
8, 300(1^/2.5 kpc) km s~^, where D is the distance (Vancura et al. 1995; Hwang & Gotthelf 
1997; Hughes 1997). X-ray spectra of Tycho's remnant further suggested that the Fe fine 
emitting gas is in general at a higher temperature and lower density than than the Si line 
emitting gas (Hwang et al. 1998); this result is also expected from the Ni bubble effect. 
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Wang & Chevalier (2001 and 2002, hereafter WCOl and WC02) studied the role of 
ejecta clumps in the evolution of SNRs. The ejecta knots associated with Tycho's remnant 
require a free expansion velocity v ~ 7, 000 km s~^ and a density contrast x ^ 100 relative 
to the surrounding ejecta, to survive crushing and remain distinct features in the remnant 
in the present epoch (WCOl). The remarkable protrusions (Aschenbach et al. 1995) on the 
periphery of the Vela SNR are also likely to be caused by ejecta clumps, and we estimated 
that X ~ 1000 ^iid ~ 3,000 are needed to create the structure (WC02). My previous 
work (Wang 2005, hereafter W05) further examined the process of ^^Ni radioactive heating 
as a plausible mechanism for forming such ejecta clumps, taking into account the effect of 
radiative diffusion. In the case of SNe la, the created clump properties seem compatible 
with those needed by Tycho's knots. This result nevertheless contrasts with the case of 
core-collapse SNe like Vela's bullets, whose required high compression is not anticipated in 
the simplest Ni bubble scenario. 

Evidence may exist that for SNe la, ejecta clumping occurs around one large central Ni 
bubble. In radio observations of Tycho's remnant, Reynoso et al. (1997) and Velazquez et 
al. (1998) found a wide and regularly-spaced emission in the NW edge of the remnant, which 
extends very close to the blast wave. The large radius of the structure is not attributable to 
the intershock instabilities that grow from small perturbations, but is effectively explained by 
the expansion of a spherical shell of clumpy ejecta into the intershock region of the remnant 
(WCOl). In fact, compared with various extensively observed Galactic Type II SNRs which 
in general exhibit quite complex abundance structures (such as in Hughes et al. 2000 and 
Hwang, et al. 2000 for Cas A; and Park et al. 2002 for G292.0+1.8), both observations and 
modeling of X-ray emissions of Type la SNRs revealed a simple and stratified composition 
(Badenes et al. 2003, 2005; Decourchelle et al. 2001 for Tycho; and Lewis et al. 2003 for 
N103B). The stratification of ejecta of SNe la seems to be a natural consequence of the lack 
of neutrino-driven convection (Kifonidis et al. 2000) and pulsar winds (Blondin et al. 2001b) 
in them, which two mechanisms can mix elements and enhance compression. Conversely, 
the disturbance that is expected in core-collapse SNe may greatly complicate the evolution 
of the Ni bubble on a large scale. 

Dwarkadas & Chevalier (1998, hereafter DC98) reviewed the solutions for the density 
structure of SNe la obtained by various 1-D explosion models; they found that an exponential 
profile, which effectively evolves from a steep power law profile to a fiatter one, is the best 
overall approximation to the density profile that is obtained by summing over individual 
chemical elements. However, in W05, a uniform density distribution (comparable to the 
inner flat component of an n = 8 power-law ejecta density model with the density p oc r"" 
in the outer parts) was assumed, so the results may not be suitable for SNe la. Furthermore, 
the W05 models of SNe la show that the propagation of the Ni bubble is only slightly faster 
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than the free motion of the ejecta, suggesting that the case of SNe la is more susceptible 
than the case of core-collapse SNe to the ejecta substrate structure. However, how the actual 
exponential decline in density affects the solution is unclear. 

The aim of this work is to determine the properties of the Ni bubble-driven shell as 
ejecta clumps associated with SNe la, using a realistic ejecta structure. A simple, adia- 
batic and spherically symmetric scenario is first considered to determine how an evolving 
exponential density gradient influences the predicted clump properties. Radiative-transport 
radiation hydrodynamics (RHD) arc then adopted to include the effect of the diffusion of 
radiative energy across the Ni bubble structure to show how the radiative solution differs. 
The explosion parameters used in our calculations arc similar to those of the 1-D explosion 
models of Hoflich & Khokhlov (1996, hereafter HK96) and Hofiich et al. (1998, hereafter 
HWK98). We note that the radiative transport process due the ^^Ni radioactive decay has 
not before been incorporated in the hydrodynamic (HD) phase of the explosion models, and 
such a proper treatment is particularly difficult for numerical convergence in the case of SNe 
la . § 2 elucidates our computational setup and methods. § 3 presents the evolutionary 
properties of the Ni bubble shell. § 4 draws on the self-similar solution for a pulsar bubble, 
to provide insight into the shell structure. § 5 compares the inferred ejecta clump properties 
of various exponential models. § 6 draws conclusions. 



2. DENSITY STRUCTURE IN TYPE la SUPERNOVAE 

SNe la are widely regarded as being thermonuclear explosions of carbon-oxygen (C-0) 
white dwarfs (WDs). Soon after the explosion, the ejecta freely expand so that each gas 
clement moves with a constant velocity v = r/t and its density drops to t~^, as in a spherical 
expansion. DC98 showed that, in this phase, the density distribution of valid SN la explosion 
models can generally be described by an exponential profile, 

PSN = Aex.p(-v/ve) r^, (1) 

where A is a constant and Ve is another constant called the velocity scale height, which are 
determined by the total mass M and kinetic energy E of the ejecta. The two constants 
can be obtained by integrating the mass density and the kinetic energy density over space: 
M = SnAvl and E = ASnAvl, or 

Ve = (£;/6M)^/2 = 2.44 X 10^ M~]/^ Egf cm s'^ (2) 



(3) 
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where M1.4 is the explosion mass in terms of the Chandrasekhar mass, IAMq, and £'51 
is the explosion energy in units of 10^^ ergs. A comparison with the power-law model 
p oc v~'^t^^ (Chevalier & Liang, 1989) yield the approximate power index of the density, 
n = —dlnp/dlnr = vjv^ = rjv^t. Hence, models with higher velocity scale (increasing with 
the E/M ratio) yield a flatter density distribution at a given flow velocity, and the gradient 
of the density increases with radius and decreases over time. 

The creation of an exponential proflle for SNe la is probably related to the fact that 
the explosion energy is steadily released behind a burning front, unlike in the case of pure 
shock acceleration in core-collapse SNe. Three progenitor scenarios have been proposed 
to distinguish SNe la (HK96; HWK98; for summary see Hoflich et al. 2003). The most 
favored is the explosion of a near-Chandrasekhar mass C-0 WD that has accreted mass 
through Roche-lobe overflow from an evolved companion star, triggered by compressional 
heating near the WD center. In this scenario, ignition is initiated in the central high-density 
region of the WD. In the deflagration model, the burning front propagates subsonically 
through the white dwarf; a popular example is the W7 model of Nomoto et al. (1984). In 
the other class of models, which includes delayed detonation (DD) (Khokhlov 1991) and 
pulsating delayed detonation (PDD), the initial subsonic burning front gradually makes a 
transition to a supersonic burning front (detonation). In the DD models, the detonation front 
erases the chemical structure that is left behind by the deflagration and creates a layered 
chemical structure, as observed. With a particularly narrow range of energy variation, the 
DD models reproduce the optical and infrared light curves of typical SNe la quite accurately. 
In current 3-D pure deflagration models, however, a signiflcant fraction of the C/0 WD 
remains unburned and mixes with the burned material at all velocities, in disagreement with 
observations. 

In the second scenario, or DET2env models, combustion occurs in a rotating low-density 
WD that is surrounded by an extended C-0 envelope. The conflguration forms during the 
merger of two low-mass binary WDs which lose angular momentum to gravitational radiation. 
The combined mass of the system may exceed the Chandrasekhar limit, but neither WD is 
above the limit. The third scenario is a double detonation of a sub-Chandrasckhar mass C-0 
WD, or the He detonation (HeD). The explosion starts with the detonation of a helium layer 
that has accreted from a close companion, which subsequently triggers the detonation of 
the inner C-0 core. The exponential model is a good simpliflcation of the density structure 
obtained from various 1-D explosion models, as presented in Fig.l of DC98. Both the 
W7 model and the delay detonation model DD200c yield a perfectly exponential decline 
in density distribution, and they predict similar properties (HKW98). The PDD and the 
HeD (PDD3 and PDDlc; HeDlO and HeD6) also exhibit more of an exponential than a 
power-law density distribution, although the gradient appears to drop in the innermost 
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regions. DC98 commented that the flattening may have been exaggerated by the numerical 
effect of inadequate zoning near the center. In the PDD and some of the HeD models, the 
explosion energy is below 10^^ ergs. As low-energy models tend to be associated with a high 
density gradient, the density evolution of lower-energy models is more sensitive to numerical 
resolution. The merger scenario is believed to be responsible for only a small fraction of 
the SNe la population because the predicted large amount of unburned C/0 at the outer 
layers is inconsistent with the IR spectra. Though not considered in DC98, this scenario may 
likewise bear an exponential density decrease, since a pure thermonuclear runaway reaction 
(as opposed to a subsonic burning) would incinerate almost all of the WD to form ^^Ni, 
which process would be inconsistent with the observations of SNe la that only ~ 0.6Mq of 
s^Ni is generated (Hoflich 2003). 

The chemical structure of SNe la varies with the explosion mechanism (HK96; DC98; 
HWK98; Hoflich 2003). In the accretion scenario, thermonuclear runaway takes place near 
the center, such that neutron-rich isotopes such as Fe, rather than ^^Ni, are synthesized 
in the core. However, for a model such as DD200c (HWK98), the central void of ^^Ni 
amounts to only ~ 10% of the total synthesized mass. The W7 model, whose density 
structure can be used as a prototype for an exponential model, shows that the amount of Si 
substantially exceeds that of Fe in the denser parts of the shocked interaction region, while 
the Fe dominates close to the reverse shock, for the age and ambient density of Tycho's 
remnant (DC98). A similar composition is produced in delayed detonation or pulsating 
delayed detonation models for Tycho's remnant, except that the densest parts of the shocked 
ejecta appear to consist mainly of C and/or O. In the merger scenario, by contrast, burning 
occurs in the low-density regions, causing little neutronization; therefore, ^^Ni is synthesized 
in the innermost region. In the HcD scenario, two separate layers of ^^Ni are generated. The 
outer one is formed by burning in the He envelope above the velocity range 11,000-14,000 
km s~^. The mass of the inner ^^Ni can be as little as that of the outer ^^Ni, or it can be 
considerable and hence represent a bright SN la such as, for example, model HeDlO of HK96. 
In the models HeD6 and HeDlO, the outermost dense regions of the shocked gas comprise 
mainly Fe, with O dominating in a dense shell inside Fe. The HeD scenario was advocated in 
view of the large population of low-mass stellar progenitors. However, Hoflich & Khokhlov 
(1997) maintained that its spectra are generally too blue to reproduce the observed hght 
curves and late time spectra. 

Subluminous models giving acceptable fits to the observed SNe la; for example, those 
in HK96 using a sample of 26 SNe la tend to eject similar amounts of Fe along with the 
^^Ni. In such models, the initial mass of ^^Ni ranges from 0.1 to 0.18 Mq, while in normal 
models, it ranges from 0.49 to 0.83 Mq. The properties of the Ni bubble expansion are 
dominated by the initial velocity distribution of ^^Ni and the effect of thermal pressure is 
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negligible (§ 3). Therefore, including a realistic ejecta elemental distribution in our models 
equivalently reduces the initial density contrast of ^^Ni, which is assumed to be distributed 
over the whole SN interior, with a mass fraction equal to uj everywhere within the bubble, 
such that UJ < 1. The range 0.5 < a; < 1 is presumed. Such a simple parameterization 
of the ^^Ni distribution may sufficiently reproduce results based on a more detailed ^^Ni 
distribution (§ 3.4). 

Figure 1 of DC98 shows that the explosion models significantly diverge from the expo- 
nential models as the velocity exceeds > 10, 000 km s~^, which transition presumably marks 
the burning interface between the core and the unburncd envelope. Yet, inside the interface, 
the drop in density appears very smooth, and the results suggest that the Ni bubble is not 
capable of expanding close to this region. 

To describe the evolution of an SNR in an exponential model, a set of scaling parameters, 
R', V and T', are used following DC98 and WOl: 



pa 2.19 

r' = ^;« 248 B-'^(|-)^'\-/3yr. (6) 

where is the ambient ISM density of the SNR and Ham = Pam/(2.34 x 10~^^) gm cm~^, 
suitable for an ISM with an H/Hc ratio of 10/1 by number. The dimensional variables r, 
V and t can be expressed in terms of the non-dimensional quantities r' = r/R', t' = t/T' 
and v' — v/V. Since the deceleration of the SNR is caused by the ambient ISM, a higher 
surrounding density is responsible for a later evolutionary phase of the remnant (larger T") . 
Once the non-dimensional solutions are obtained, they can be re-scaled to the dimensional 
solutions for a different set of M, E and pam- Thus, one evolutionary sequence in the 
non-dimensional variables represents virtually all possible dimensional solutions. 

Figure 1 plots the evolution of the characteristic velocities of an SNR that arises from 
the interaction of an exponential ejecta with a constant ambient medium, corresponding to 
Fig. 1 of WCOl. As the ejecta density profile becomes fiatter, the reverse shock, initially 
moving outward in the stellar frame, begins to turn inward at t' ~ 2.5. It reaches the 
stellar center at t' ^ 8. Unlike in the power-law model in which the evolution is self-similar 
(Chevalier et al. 1992), the exponential model gives rise to different deceleration rates as the 
remnant evolves. Nevertheless, the fiow velocities immediately below the reverse shock front 
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do not substantially deviate from those in the power-law model (Fig. 12 of W05) within the 
epoch t' ~ 0.1 — 2', when the clump-remnant interaction is likely to begin. 

3. NICKEL BUBBLE SHELL STRUCTURE AND EVOLUTION 

3.1. Methods 

The computational methods used in this paper are similar to those adopted in the 
author's previous study W05. HD and RHD simulations were performed based on the two- 
dimensional finite-difference code ZEUS2D (Stone & Norman 1992; Stone et al. 1992). 
Heating due to a central core of ^^Ni that has a density contrast ui = 1 relative to the 
ambient ejecta was considered. The ejecta were initially freely expanding with an exponential 
density profile that is characterized by an explosion mass M and energy E. The radioactive 
energy was deposited continually to the gas thermal energy (HD case) within the bubble 
as the bubble-shell interface was tracked in a uniformly-expanding grid in spherical polar 
coordinates, starting at to = 100 s after the explosion. The inner boundary of the grid had 
a simple refiecting condition. At the outer boundary, a non-zero gradient outfiow condition 
was applied to preserve the ambient exponential profile and eliminate spurious shocks that 
were caused by the grid expansion. Basko (1994) and W05 described the radioactive power 
input. Since the estimated mean life of ^^Ni is larger than that suggested in Firestone (1996), 
< 2% less radioactive energy than indicated by Firestone was accumulated in the radiation- 
streaming epoch ~ 10^ s. A ^^Ni abundance of 0.5Mq yields an effective energy input of 
~ 6 X 10^^ erg, or 6% of the initial kinetic energy. Therefore, the increase in the velocity 
from local free expansion is typically only around < 3%. 

The initial ejecta velocity at the bubble edge was determined from the initial ^^Ni density 
contrast and age: C/q = Ro/to ~ u;~^/^. The background thermal pressure was distributed 
according to p = np^, where 7 = 4/3 in the HD case, and k — 6.1x 10^^ (cgs units) (derived 
from the change in entropy during nuclear burning of C-0 to ^^Ni, W05). The late-time 
radioactive power input prior to radiation streaming dominated the pressure in the bubble. 
In the RHD simulations, the radioactive power was deposited into the radiative energy and 
7 = 5/3 for matter. In supernovae, since the dynamic time required for the ejecta to become 
freely-expanding is much longer than that required for the radiation field to be thermalized 
with electrons and for collisions between electrons and ions to occur, a black-body radiation 
field is well established (Shigeyama & Nomoto 1990). To ensure proper evolution of the 
flow toward the optically-thin regime, full radiative-transport RHD is invoked, under the 
assumptions of local thermal equilibrium (LTE) and gray opacity. The opacity of the ejecta 
is dominated by electron scattering; absorption is neghgible (Shigeyama & Nomoto 1990). A 
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Thomson scattering cross section of 0.2 cm^/gm was used to calculate the scattering opacity, 
which is suitable for a gas of completely ionized heavy elements. The time for radiative 
energy to diffuse out of the bubble depends on the extent of electron scattering and the 
column density of the shell. Notably, full-power radioactive heating can not be achieved in 
the radiative case because radiative diffusion limits shell expansion. 

3.2. Hydrodynamic Simulations 

Figure 2 and 3 show the ejccta structure and the evolution of the dynamic properties 
of the Ni bubble shell in the adiabatic case, with parameters M = 1.4 Mq, E = 10^^ ergs, 
M^i = 0.5 Mq and c<j = 1, in what is referred to herein as the standard model. At the 
commencement of the simulation to = 100 s, the bubble has the radius Rq = 5.18 x 10^° cm, 
velocity Uq = 5.18 x 10^ cm s^^ and ambient density 0.92 g cm~'^ at its edge. The initial 
background thermal energy density is 2.77 x 10^^ erg cm~^ in the center. By this stage, a 
total radioactive energy of 4.85 x 10^^ ergs has been deposited since the explosion, yielding 
an averaged radioactive energy density of 6.89 x 10^^ erg cm~^ in the bubble. 

The inflation of the bubble induces a strong forward shock behind which the ambient 
gas is compressed into a dense shell (Fig. 2). A notable feature of the structure is the drop 
in density in the ejecta substrate. The density contrast across the shock front is seven, as 
for a radiation dominated (7 = 4/3) strong shock. The inner edge of the shell is a contact 
discontinuity where the gas has been shocked and cooled for the longest time and so has the 
highest density. To describe the shell acceleration, the shock radius is approximated as a 
power-law function of time, Rgh oc t"-, where a is the expansion rate at the shock front, which 
is equivalent to the velocity contrast between the shock and the ambient freely-expanding 
ejecta, {dRsh/dt)/{Rsh/t). The expansion rate rises to the maximum a < 1.05 around 10^ s 
and then falls to 1.0. That is, the shock front first accelerates and then freezes in the local 
free expansion of the fiow (Figs. 3(a) and (b)). The shell cannot be slower than the free 
expansion. 

The shell is very thin; at 10^ s it reaches a maximum thickness ratio of /3 = hgh/Rsh ^ 
0.004, where hgh is the thickness of the shell (Fig. 3(c)). A linear relationship holds be- 
tween the thickness ratio and the expansion rate, /3 ~ (a — 1)/10, because the expan- 
sion rate is low. The Mach number of the shock, M — {dRgh/dt — Rgh/t) (po/lPo)^^^ — 
(a — l)Rsh/t (po/7Po)^^^, is Al~ 87 at 10^ s, where po and po are the preshock density and 
pressure, respectively. The expansion proceeds until pressure equilibrium is reached, or in 
the radiative case until the shell becomes transparent to the 7-rays. The shell's compression 
ratio to the preshock ejecta is x > 10 (Fig. 3(d)), and is limited by numerical resolution. The 
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compression in the bubble interior does not drop with time, as it does in the uniform-density 
model; rather it increases to ~ 0.9 (Fig. 3(d)). 

At ~ 10^ s, the shell has a thickness ratio (3 > 0.002, a velocity Vf ~ 6,000 km s~^, 
and a swept-up mass Mg ~ 0.09 Mq. At this time, the surface density of the shell begins to 
drop below the mean free path for the ~ IMeV 7-ray photons, ~ 33 g/cm^, determined by 
Monte Carlo simulations of X-ray transfer in SNe la (corresponding to an effective opacity 
of 0.03 cm^ g^^; Sutherland & Wheeler 1984, Shigeyama & Nomoto 1990). The gas is then 
transparent to the radiation. Although 75% of the total radioactive energy is yet to be 
released (Fig.l of W05), further deposition of the radioactive energy does not affect the 
dynamics of the flow, and the flow evolves to a free expansion. 



3.3. Variation of Initial Parameters - Adiabatic Caises 

The effect of explosion parameters and other initial conditions on the solutions is now 
examined. A wide range of explosion parameters, such as those used in the 1-D explosion 
models of HK96 and HWK98, have been studied, under both luminous and underluminous 
conditions (Tab. 1). The clump strength is conservatively evaluated by approximating the 
shell's frozen-in velocity by the flow velocity of the contact discontinuity when the shell's 
surface density starts to drop below ~ 33 cm^ g~^. 

A comparison with the uniform-density model (model 'flat') indicates that in the stan- 
dard exponential model (model 'std'), the initial amount of ^^Ni was distributed at lower 
velocities, such that radioactive pressure within the bubble is larger. The shell reaches a 
higher expansion rate, and the density contrast and thickness ratio are therefore greater. 
Nevertheless, the expansion velocity and swept-up mass of the shell are lower (Figs. 3(e) 
and (f)). In the adiabatic uniform model, the frozen-in velocity is Vf ~ 7,500 km s~^ and 
the swept-up mass is Mg ~ 0.12 Mq, while in the exponential model, Vf ~ 5,800 km s^^ 
and Ms ~ 0.09 Mq. The difference between the velocities of the radiative counterparts of 
the above two models should not exceed the difference between their maximum velocities, 
~ 8, 800 km s-^ and ~ 7, 500 km s-\ 

For a lower initial ^^Ni density contrast uj = 0.5, the frozen-in age is slightly advanced, 
with a lower expansion rate, a lesser shell mass M^, ~ 0.03 Mq and a higher frozen-in velocity 
Vf ~ 9, 300 km s~^. The drop in the expansion rate and the swept-up mass follow from the 
lower radioactive pressure within the bubble. The frozen-in velocity is higher because the 
^^Ni initially extends to higher velocities. In the exponential model, the radioactive pressure 
dominated the thermal pressure since the initial stage of heating; the thermal pressure thus 
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has a negligible impact. In the uniform model, reducing u results in a greater swept- up mass 
(Fig. 3(f)), because our assumed thermal pressure equilibrium across the bubble causes the 
initial thermal energy in the bubble to be higher when uj is lower. Notably, the c<j = 0.5 
standard model yields a higher shell velocity than its uniform-density counterpart, because 
of the initially higher velocities of ^^Ni. 

Comparisons of various models show that models with a higher velocity scale height Vg, 
or a flatter density distribution, tend to yield a higher shell velocity, a lower expansion rate 
and a lower swept-up mass. For example, the HeDlO model (in which the density decline is 
the least among all of the models) yields a shell velocity that extends to Vf < 13, 000 km s~^ 
and a swept-up mass of < O.O2M0, whereas the PDDlc model (which has the greatest 
density gradient) gives Vf ~ 2, 500 km s^^ (or ~ 3, 500 km s"^ at maximum) and Mg ~ 
O.I6M0. The shell velocity is correlated with the initial velocity distribution of ^^Ni, or the 
velocity scale height, because the shell expansion rate barely exceeds the free expansion of the 
local ejecta. The expansion rate is generally expected to be under a < 1.07. The momentum 
of freely expanding ejecta that were swept up into the shell, obtained by integrating the 
momentum density over the velocity space from the initial fo final shell location, is found 
similar to the final momentum imparted to the shell. Thus, only a small fraction of the final 
shell momentum results from the bubble acceleration. Figure 2 shows that the density profile 
of the bubble interior deviates from an exponential model mostly in the regions that are close 
to the contact discontinuity, also because the expansion is weak and can be contrasted with 
that in the case of core-collapse SNe (Fig. 2 of W05). Figure 3(h) plots the time evolution 
of the approximate power index, n = vlv^^ at the contact discontinuity. The frozen-in shell 
velocity does not exceed about three times the velocity scale height of the model - except 
for the luminous HeDlO and low-o; models. 

The swept-up mass of the shell is not very sensitive to the initial abundance of ^^Ni. 
In a Chandrasekhar-mass, low-energy model using £"51 = 0.7 and Mjvi = 0.49, for example, 
the swept-up mass is comparable to that of a similar model with half the amount of ^^Ni. It 
also is approximately the same across several models with intermediate explosion conditions, 
with a typical Mg of ~ O.IM©. 

Solutions to the underluminous models are not sensitive to variations of the initial 
parameters including the initial ^^Ni density contrast, because their steep density distribution 
and low ^^Ni abundance allow little corresponding variations in velocity and radioactive 
pressure. A unique feature of the exponential models except for PDDlc is that the density 
contrast at the bubble interior increases over time, because of continual flattening of the 
exponential density distribution and a relatively low expansion rate. 

In the HeD scenario, the outer layer of ^^Ni is present at high velocities of above 11,000- 
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14,000 km (HK96). An off-center fieating scenario was considered using tfie HeD6 model 
(model 'HcD6*'), in which the radioactive power from a half of the outer ^^Ni mass was 
gradually deposited in the space between the velocities 12, 000 km s~^ and < 13, 400 km s~^, 
while the inner grid boundary continues to move at a fixed velocity of 12,000 km s^^. A 
dense shell as in the central-heating model is swept up, but it exhibits very weak acceleration 
relative to the local fast ejecta. The onset of the optically-thin stage thus advances to ~ 10^ s, 
leading to a 'low' frozen-in velocity at < 14, 000 km s~^. 

A fiattcr distribution of the thermal pressure, p ^ p, was assumed with a constant 
temperature, kT ~ 0.6 MeV (derived from the nuclear burning process, W05). As the 
central background pressure is enhanced by a factor of ~ 500, the radioactive pressure still 
dominates the thermal pressure, so the evolution of the Ni bubble remains largely unchanged. 
To account for the central ^^Ni void as in the accretion scenario, ^^Ni in model DD200c was 
redistributed to higher velocities, corresponding to a void mass of 1/10 of the synthesized 
mass. The result is indistinguishable from that of the central-heating model. Lastly, the 
simulated solutions that began at the age to = 1000 s rapidly converged with those that 
began at 100 s. Thus, as long as to ^ 1000 s, the final solutions do not depend on the 
commencing ages of the simulations. 



3.4. ^^Ni Mass Fraction in the W7 Model 

In our u = 1 exponential W7 model, the shell velocities, 7558 km s~^, is actually 
lower than the minimum velocity ~ 10, 000 km s~^ of the Si-dominant layer, obtained using 
the W7 model of Nomoto et al. (1984). In their original calculations, the velocity of the 
radioactive Ni was as high as 8300 km/s, and its fractional mass fell to 0.5 at 8800 km/s. To 
justify that a simple parametrization of the initial ^®Ni distribution in terms of 0.5 < a; < 1 
approximates the results that consider a reahstic ejecta elemental distribution, the ^^Ni mass 
fraction given by the model of Nomoto et al. {Mjsfi — 0.58Mq) was incorporated into the 
exponential model. As the ^^Ni mass fraction falls steeply above 8,800 km s~^, the initial 
^^Ni distribution was cut off at this velocity but the same radioactive power, proportional 
to the mass density within it, was deposited. Spatial variations in the background thermal 
pressure caused by the nonuniform elemental distribution were neglected. This modified 
model gives a similar initial outbound speed to that obtained in our uj = 0.7 W7 model 
(Fig. 4). (The integrated density over all elements should follow an exponential profile as 
illustrated in DC98.) In both models, the final shell velocity reaches ~ 9000 km s^^. In the 
modified model the swept-up mass is lower but does not differ significantly from that in the 
CO — 0.7 case. 
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The reason why the use of a; = 0.7 adequately reproduces the result based on a more 
detailed ^^Ni distribution in the W7 model is that it places the outer boundary of ^^Ni 
in the 8300-8800 velocity range, beyond which whose mass fraction drops down rapidly. 
However, in view of a finite sound speed propagation within the bubble and the presence of 
nonradioactive material in the ejecta, this result strongly depends on the assumed value of 
uj and is merely an approximation. On the other hand, the clump speeds in our uj = 1 cases 
stand only as the lower limits in more realistic Ni bubble heating scenarios. 

3.5. Radiation-Hydro dynamic Solutions 

Here we present full-radiative transport RHD simulations to illustrate the effect of the 
diffusion of radiative energy across the bubble structure on the shell properties. In Figs. 2, 
4 and 5, the density structures of the radiative models are plotted over their adiabatic 
counterparts. 

The shock has a widely extending radiative precursor, lowering the overall compression 
of the shell. When the radiative precursor is considered, the shell mass reaches ~ 90% of the 
adiabatic value (0.1 Mq). Excluding the radiative precursor reduces the mass to > O.O4M0 
(radiative case) in the standard model. The drop is typically < 50%. 

Since the acceleration of the shell is weak, the shell formation processes can be regarded 
as perturbations to the initial free expansion, such that the final radius and velocity of the 
shell are only slightly reduced. The final shell velocity differs from the adiabatic solution by 
< 3% (Fig. 6). As the final shell velocity may not be an adequate measure of the radiative 
effects, the difference between the final shell velocity and the velocity of the initial free 
expansion at the bubble's edge is used to assess the importance of radiation diffusion effects. 
A ~ 25% difference is found between the adiabatic and radiative solutions. However, in the 
cj = 0.7 W7 model, it is only 12%. Hence, the radiative effects are expected to be mild in a 
more realistic scenario. 

The mass and number of clumps should drop as the acceleration and compression of the 
shell are decreased. Since the RHD simulations are affected by inadequate spatial resolution, 
the shell's density contrast is crudely estimated, based on the assumption that that the shell 
mass is 50% less and the shell thickness is double the size of that in the adiabatic solutions; 
a compression ratio of ~ 25% of the adiabatic value, or x ~ 25, is expected. Nonetheless, 
since ejecta clumps (shell fragments) most likely arise from the densest part of the CD, the 
true compression ratio may be higher. 
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4. SELF-SIMILAR APPROXIMATIONS 

To provide insights into the exponential model, this section discusses the self-similar 
solutions of Chevalier (1984) and Jun (1998) for a pulsar bubble that arises from a pulsar 
wind that interacts with supernova ejecta. The analogy to the problem of interest follows 
from the fact that the shock is driven by the motion of CD, as in the process by which a 
piston shoves up mass in a gas chamber. Therefore, the properties of the forward shocked 
layer (shell) can be determined by the motion of the piston (CD) and the density of the 
gas downstream (outer ambient ejecta); details of how the inner gas is heated to power the 
piston need not be known to apply self-similar solutions. 

The Chevaher (1984) and Jun (1998) solutions considered a pulsar wind with a time- 
varying luminosity L oc t"^ and supernova ejecta with a constant power index n associated 
with the denstiy decline. A simple dimensional analysis reveals the self-similar expansion rate 
of the bubble to be a^^ = (6 — n — /)/(5 — n). Accordingly, the rate of expansion increases 
with the density gradient. The self-similarity exists only for 0^,.^ > 1, or I < 1, since the shock 
cannot be slowed down in the ejecta to be slower than the free expansion, and for ^ ^ 3, 
since the mass must be finite. With the given adiabatic index 7 = 4/3, self- similar solutions 
can be described by two parameters: ass and n. Since the exponential density gradient 
increases with radius and flattens over time, n is evaluated at the contact discontinuity near 
the frozen-in age of 10^ s (Fig. 3(h)). In our standard case, n — 2.4. 

In the evolution, the shock is expected to accelerate (a > 1) to the self-similar regime, 
and so the shell thickness rises continuously. If a constant luminosity I — and n — 2A 
powers the expansion, then Ugs — 1.38. In our hydrodynamic solution (simulated solution 
of the standard model), the shell thickness ratio reaches a maximum at 5 x 10^ s; this is 
correlated with an increase in the accumulated radioactive energy at 5 x 10^ s (Fig. 4 and 9 
in W05). The turnover after 10^ s follows from the fact that / > 1 (where / is estimated by 
taking the time-derivative of the accumulated radioactive energy, / ~ —dlnq{t)/dlnt, where 
q{t) is as given in W05), when the self-similarity breaks down; it is not related to the mass 
of ^^Ni or the gradient of the ejecta. The decrease in the shell thickness during the earliest 
evolution is due to the negligibility of the medium motion during a sudden acceleration (Jun 
1998), which can be regarded as a numerical artifact that depends on the starting age of the 
simulation. 

The hydrodynamic solution has a weaker expansion than the self-similar solution for 
a fixed density gradient n — 2.4, but gradually approaches the self-similar regime. In the 
self-similar solution, ass — 1-07 at 1 x 10^ sec (for I ~ 0.8, which is the luminosity index at 
1 X 10^ sec), and ags — 1.0406 at 2 x 10^ s (for I — 0.896), while in the hydrodynamic solution, 
a — 1.04 and 1.0444 at the respective ages. The self-similar solution thus approximates to the 
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hydrodynamic solution near the frozen- in stage (Fig. 7). The self-similar density is infinite at 
the contact discontinuity. Chevalier (1984) gave the asymptotic density profile close to the 
contact discontinuity in the uniform model: p oc {j'jRcD — 1)"^, where h = 2ass/{7ass — 5). 
This result can readily be applied to a power-law model. Because of the sharp increase in 
density toward the CD, the coarseness of the grid obviously limits the highest computed 
density in the grid domain. 

The self-similar solutions to the thickness ratio f3 and average density contrast Xavg of 
the shell are derived as a function of the self-similar expansion rate a^g (Fig. 8). Both /? and 
Xavg increase with the density gradient. Xavg reaches a maximum (turnover) with respect to 
Gss in cases of n < 3. For example, for n = 2.4, the density contrast turns over at ass = 1-04 
and the shell thickness ratio is about 0.004. /3 is under < 0.023 when the largest expansion 
rate allowed for an expanding medium, Uss — {6 — n)/{5 — n), is achieved. (The Chevalier 
and Jun solutions used gamma=5/3 and thus yielded a larger thickness ratio, about 0.02 for 
n — and a = 1.2.) In the n = 2.4 case, /3 ranges from ~ 0.0055 to ~ 0.006 for ass ^ 1-06, 
and Xavg ^ 75, depending only mildly on the expansion rate between 1.02 and 1.08. Xavg 
may be enhanced by a factor of ~ 50% over that in the n = case. Figure 7 displays a 
larger shell thickness in the hydrodynamic solution than in the self-similar solution. If the 
hydrodynamic solution reaches an equally infinite density at the CD, then the true average 
density contrast of the shell may be even higher than that given by the self-similar solution. 

We caution that, the self-similar solution can also be affected by the spatial resolution 
used in integrating the governing self-similar differential equations from the shock front to 
the CD. High-resolved results are generated using the fourth and fifth -order Runge-Kutta 
methods with a small integration step size until the density contrast values do not increase 
further (Fig. 8). Fluctuations in the curves are caused by the accumulation of errors in 
summing up the shell mass over the acquired density distribution. 

We have shown that the hydrodynamic solution of an exponential model can be approx- 
imated by the self-similar solution of a power-low model. Nonetheless, such a description is 
merely semi-quantitative for the problem under consideration, as it has no true self-similar 
solution. 



5. NICKEL BUBBLE SHELL AS EJECTA CLUMPS 

This section examines whether the ejecta clumps that originate in the Ni bubble shell 
can initiate the clump-remnant interaction sufficiently early in time and further survive 
crushing in a remnant, such as Tycho's knots, under various explosion conditions (Tab. 1). 
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As a shock moves past a cloud, the incident shock generates a transmitted 'cloud shock' that 
moves into the cloud and crushes it. When the cloud shock exits the cloud, a rarefaction 
wave moves back into the cloud and causes lateral expansion. The combined instabilities 
destroy the cloud on a timescale of several times the cloud-crushing time, t^c = ''"cX^^^/'^Oi 
where Tc is the cloud size and vq is the initial shock velocity (Klein et al. 1994). When a 
clump interacts with a remnant, the clump survival depends on three parameters: (1) the 
initial density contrast of the clump, relative to the surrounding ejecta; (2) the initial 
impact age of the clump with the reverse shock, Irs] this is related to the absolute ejection 
speed of the clump, which is given by the frozen-in velocity of the Ni bubble shell, as well as 
the true age of the remnant, and (3) the initial size of the clump. Of the three parameters, 
the clump speed and initial impact age with the remnant are barely limited by resolution and 
these properties are well-defined for comparing observations. As they do not significantly 
diverge in the case of radiative diffusion, the following discussions are primarily based on the 
adiabatic solutions. 



5.1. Tycho's SNR as the Standard Model 

Tycho's remnant is used as a prototype for the clump-remnant interaction. Table 1 lists 
the speed of the clump (frozen-in velocity of the shell) for each model under consideration. 
The initiating age of a clump-RS impact when the clump speed matches the fiow velocity 
immediately below the RS is derived using the velocity evolution of a remnant in Fig. 1, The 
dimensional values of the ages are noted specifically for two ISM density values, Uam = 0.5 
and 1. The final column of Table 1 presents compares of the normalized age of Tycho's 
remnant (whose real age is 435 yrs) with the RS-impact; the comparison yields favorable 
explosion conditions for Tycho's supernova. 

Since the remnant of a Type la SN continuously decelerates with time, its explosion 
parameter and ambient density can be inferred from the observed deceleration rate of the 
remnant, which parameter is independent of distance, given its dimensional age. For a 
Chandrasekhar-mass, E^i — 1 explosion, DC98 noted that the observed deceleration rate 
5 ~ 0.5 (DC98 and references therein; Hughes 2000) of Tycho's remnant depends on an ISM 
density in the range of Uam = 0.6 — 1.1 cm^^. In this case, Tycho's dimensionless age is 
between t' ~ 1.39 (for = 0.5 cm^"^) and ~ 1.75 (for Uam = 1-0 cm^^). For Uam = 1, 
the reverse and forward shocks have radii of 2.08 pc (r' = 0.98) and 3.09 pc (r' = 1.45), 
and decelerate with S = 0.15 and 6 = 0.47, respectively; the result is consistent with the 
observations. However, since the deceleration of the forward shock evolves very slowly near 
Tycho's present age, t' ~ 2.0, when 6 ~ 0.5 (inset of Fig. 1), the explosion parameters can 
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be rather unconstrained. 



5.2. Clump Speed and RS-Impact Age - Indication of Explosion Conditions 

As our standard cu = 1 model yields a lower clump speed than the uniform model (for 
which t'j^g ~ 0.9), the onset of a clump-remnnat interaction is delayed to the age t'j^g ~ 1.3, 
or ~ 330 yrs, for Uam — 1 (see Fig. 1). The clump- RS impact still precedes Tycho's present 
epoch, but the clump-FS impact does not, and is anticipated to occur at t'pg ~ 2.2, or 
~ 545 yrs (Fig. 9). The clump speed is ~ 30% short of the velocity of Tycho's knots, 
V ~ 8,300 km s~^, estimated at a distance of 2.5 kpc. Nonetheless, measurements on 
Tycho's distances range from 1.5 to 4.5 kpc (see Schaefer 1996, Ruiz-Lapuente 2004 and 
references therein), so the observed velocity of the knots can also be quite uncertain. 

Explosion conditions similar to those in the DD200c and W7 models result in a timely 
clump-remnant interaction as well as a compatible clump speed. The luminous HeDlO model 
gives a high clump speed but is not favored as its clump-remnant interaction {t'^g = 0.56 
and t'pg ~ 1.0) may be initiated too early to be sustained to Tycho's current state. The 
HeDlO model is expected to yield a weaker clump strength than the typical models because 
of its lower shell expansion rate and density contrast. 

The underluminous PDDlc model yields a low expansion velocity ~ 2500 km s~^ where 
the exponential model becomes invalid (Fig. 1 of DC98). If a constant cjccta density of 
10~^^ g cm~^ at 10^ s or 1 g cm~^ at to — 100 s is assumed in this velocity space, then the 
initial bubble velocity is 3600 — 4500 km s"^ for a; = 1 — 0.5, and the final shell velocity 
may be 4000 — 5000 km s~^, which is still insufficient for Tycho's knots. For models with 
a moderate clump speed, the high-mass DET2env6 model notably gives a relatively late 
clump-RS impact, Irs ^ 400 yrs, for nam = 0.5 cm~^, while the low-mass HeD6 model 
yields a timely interaction at t^s ~ 270 — 340 yrs with a similar clump speed, despite a 
low ^^Ni abundance Mjvj = O.25M0. The uncertainties in the observed clump speed and 
deceleration rate of the remnant suggest that Tycho's SN may possibly have originated in 
the HeD6-like, sub-Chandrasckhar mass and low-energy explosion scenario. 

As the density of the SN ambient medium declines (e.g., = 0.5 cm~^, see Tab. 1), the 
clump-RS impact is further delayed, and the clump undergoes a longer intershock-crossing 
process (Fig. 10). The lower density limit Uam = 0.6 cm~^ for Tycho's deceleration rate 
in DC98 could however be underestimated, as the undecelerated motion of clumps behind 
the blast wave was not considered. Given that the dimensionless time scale T' increases by 
a factor of 2.15 for a ten- fold drop in ambient density, if the minimal ambient density is 
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0.05 cm ^, then the RS impact should occur within an age of ~ 2000 yrs under all plausible 
explosion conditions of SNe la. 

In the case of an initially lower ^^Ni density contrast a;, a faster but smaller and less 
dense clump is created. Consider the standard cu — 0.5 model as a reference: the cloud- 
crushing to intershock-crossing time ratio reaches only < 2%; the shell to Si knot mass ratio 
falls to ~ 1/13, and the forward shock impact advances to t'pg ~ 1.0. The clump speed 
(approximately the speed of Si) obtained using Nomoto's W7 model may be approximated 
by our uj = 0.7 W7 exponential model. The uj = 0.5 cases are particularly applicable to 
underluminous explosions which tend to eject an amount of Fe that is comparable to that of 
Ni. However, since the clump speed in such scenarios depends less sensitively on variations 
of the initial conditions, the increase in the clump speed relative to the cu — 1 case does not 
suffice to change our conclusions. 



5.3. Clump Density and Size - Indication of Clump Robustness 

The clump is expected to be of a size that is comparable to the thickness of the Ni 
bubble shell. In Fig. 9, the fractional clump size is estimated versus the age at which the 
clump-remnant interaction is initiated, for varying thickness ratios For /? = 0.004, the 
clump initially occupies ~ 1% of the intershock shock width, or ~ 0.3% of the forward shock 
radius, upon the time of the RS impact. If the lateral expansion of the clump is neglected, 
the uniform expansion of the clump (in which the density falls as ~ and the size increases 
as ~ t) would result in a fractional clump size in terms of the forward shock radius that 
equals the initial thickness ratio of the shell, when the clump reaches the forward shock. 

After they enter the SNR intershock region, the clumps become fragmented over a period 
of several times the cloud-crushing time (WC02, equation (8)). The self-similar solution 
indicates that the shell's average density contrast is up to ~ 50% greater than that obtained 
using the uniform model. The Rayleigh- Taylor mixing at the Ni bubble shell interface is 
thus facilitated, creating denser clumps. The cloud-crushing time scale of the clump in our 
optimal adiabatic scenario x ~ 100 and (3 = 0.004 is estimated to be tec ^ 20 yrs, or ~ 0.06 
times the intershock-crossing time {tps — ^rs)- Although the time scale seems small, in the 
high-resolved simulated clump- remnant interaction (WC02), clumps fragments successively 
create protuberances on the forward shock in the remnant, are decelerated, and the forward 
shock front catches up with them; the clump crushing process is not expedited. The dynamic 
effect of small but dense clumps thus is not negligible. Notably, for a fixed clump mass, the 

—1/2 

cloud-crushing time scales inversely with the clump size, t^c oc rc . 
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Recent Chandra observations of Tycho's SNR by Warren et al. (2005) revealed the 
presence of dense gas that stretches very close to the blast wave - much closer than the 
unstable gas according to the exponential model of WCOl. Warren et al. attributed such 
a structure to the intershock unstable gas that arises from a different intcrshock density 
profile, modified by the cosmic ray accelerations of the preshock gas. However, the small 
and dense clump fragments can also contribute to the observed dense gas; their interaction 
with the remnant does not extend the CD to larger radii. Additionally, modification of the 
intershock structure should have a weak dynamic infiuence on the dense clump fragments. 



i? 05/1 r observations of Tycho's remnant have revealed two X-ray knots that protrude 
from the SE edge of the remnant in undecelerated motion. Hwang, Hughes, & Petre (1998), 
estimated a mass of O.OO2M0 and O.OOO4M0 for the Si+S and Fe knots, respectively. The 
total swept- up mass in the Ni bubble shell {Mg ~ 0.1/0.05 Mq in the adiabatic/conservative 
radiative case) is then ~ 25 — 50 times the mass of the Si knot. While the process by which 
the shell fragments are formed is yet to be demonstrated, this mass ratio could be suggestive 
of the instability mode at the bubble-shell interface, where the density and pressure exhibit 
opposite gradients and thus satisfy the Rayleigh- Taylor unstable criterion (Chevalier 1976). 
The instability growth time can be inferred from the characteristic relaxation time r of 
a perturbation, given the differences across the contact discontinuity, in the acceleration 
(effective gravity) and position, Ag and Ax: 



By the epoch of radiation streaming 10 s, the shell has been accelerated over three orders 
of magnitudes in radius and r ~ s, so the instabilities are estimated to have fully 
developed. Alternatively, the growth of the instability can be evaluated at any stage using 
the information on the shell acceleration and thickness contained in Figs. 2(b) and (c), 
given that the growth rate of an instability mode that initially disrupts the shell is related 
to the square root of the ratio of the shell acceleration to the shell thickness, or (t/r)^ = 
(da/dlnt + a(a - I))/ (3. 

To determine the the spatial structure of SNe la, Dietrich Baade measured the polar- 
izations in the Si II 635.5 nm line, which is a very common feature among young SNe la up 
to an age of 10^ s (2006, private communication). He concluded that in order to generate 
the observed polarized effect of SNe la, there must exist a few tens (> 10) of clumps or 
bubbles in the ejecta, and their number must not be too large. However, he also noted that 



5.4. Clump Mass - Indication of Ni bubble Instability 
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the asymmetry produced in multi-dimensional SNe la explosion models yields approximately 
the same view in all directions (Reinecke et al. 2002; Schmidt & Nicmcycr 2005), so that 
none of the models can account for the polarization. To the contrary, Baade's estimate of 
the number of spatial inhomogeneities of SNe la is consistent with the shcll-to-clump mass 
ratio herein. This result indicates that the polarization may originate in the Ni bubble effect, 
which is characterized by a weak expansion. 

5.5. Indication for Tycho's SN, SN 1006 and N103B 

Indications that both Tycho's SN and SN 1006 were underluminous and corresponding 
indications for SNR N103B are finally examined. 

Van den Bergh (1993) reconstructed the light curve of Tycho's SN, which he found was 
very similar to that of the underluminous SN 1991bg (Filippenko et al. 1992; Leibundgut 
et al. 1993). Photometric calibrations and estimates of distance to Tycho's SN vary widely, 
and recent interpretations of the historical records by Ruiz-Lapuente (2004) identify Tycho 
as a normal Type la supernova. For SN 1991bg, Turatto et al. (1996) estimated an explosion 
energy that was a factor of three to five below that of typical SNe, and Mazzali et al. (1997) 
suggested a WD mass of O.62M0, a kinetic energy of £'51 = 0.62 and an ^^Ni mass of 
0.07Mq to account for the extension of the observed centrally-peaked ^^Ni to 5, 000 km s~^. 
Recent SN surveys (Li et al. 2001; Benetti 2005) classified as many as 36% of the observed 
SNe la as being intrinsically peculiar spectral types: 20% were like the overluminous SN 
1991T (Mazzali et al. 1995), while the remaining 16% were similar to the underluminous SN 
199 Ibg. As a conservative approach to understanding the inferred underluminous condition, 
a Chandrasckhar-mass and low-energy model with E^i = 0.7, Mjvj = 0.25 and u = 1 — 0.5 
were examined. The assumed Ni abundance is roughly the upper mass limit of Fe (< 
0.15 Mq) in SN 1006 (Hamilton et al. 1997). In this PDD-fike scenario, the sheU expands at 
Vf ~ 3, 800 — 5, 000 km s~^; the velocity is not consistent with Tycho's SNR. The parameters 
of Mazzafi et al. (1997) were then adopted; in this case, ~ 4000 — 5000 km s~^ for 
a; = 1 — 0.5 - roughly comparable to the observed velocity of ^^Ni in SN 1991bg. 

The peak magnitude of SN 1006 has been variously estimated at values that span over 
five magnitudes. More recent examinations by Winkler et al. (2003) imply a normal bright- 
ness, but the Fe abundance of the remnant is known to be peculiarly low (Hamilton et al. 
1997), revealing the SN likewise to be underluminous. The morphology of the SN 1006 rem- 
nant is not symmetrical; its ambient density may be three times higher toward its northwest 
region while the overall observed deceleration of the remnant is 5 ~ 0.5. DC98 proposed 
a Chandrasekhar-mass model of Uam — 0.047 cm~^ for £^51 = 1, or Uam = 0.07 cm~^ for 
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E^i = 1.3 (model 'W7-2') to explain the deceleration. With uj = 0.7 and Mm = 0.59, both 
models predict a reverse shock impact before an age of 500 yrs and a forward shock impact in 
the present epoch. With cu — 0.5 and M^i = 0.16, an expansion velocity at f / ~ 5000 km s~^ 
is yielded; this value is consistent with the lower observed velocity limit of the Si in the rem- 
nant, 5, 600 — 7, 000 km s^^, estimated for a distance of 1.8 kpc. The distance measurements 
of SN 1006 range from 0.7 pc to ~ 2.5 kpc (see DC98, Winkler et al. 2003 and references 
therein). In the case of Vf ~ 5000 km and Uam ~ 0.05 cm~^, the clump-RS impact is 
expected at an age of ~ 1000 yrs. This result indicates that the clump-remnant interaction 
may have just started in SN 1006. 

The exponential model explains the commonly-observed lack of mixing in the remnants 
of SNe la. Blondin et al. (2001a) studied the interaction of numerous small ejecta Fe bubbles 
with the intershock region of an SNR. As the reverse shock front moves back more rapidly 
through the bubbles, considerable turbulence is generated, which eventually destroys the 
original stratification of the ejecta. In their simulations, a density contrast of x ~ 1/100 was 
used between the bubbles and the ambient ejecta. Although the assumption of high density 
contrast pertained specifically to core-collapse SNe because of the expected initial mixing of 
^^Ni, such a turbulent scenario is unlikely to apply for SNe la, because the bubble interior 
has a density contrast that is close to unity, and clumps are expected to be present only in 
a thin shell at the boundary of the Ni bubble. 

A good example of the Ni bubble effect is the young and compact supernova remnant 
N103B located in the Large Magellanic Cloud (Lewis et al. 2003). In this remnant, the outer 
Si and S emissions exhibit an ionization time scale (ngt) that is ~ 100 times higher than 
that of the interior, hot Fe, while the overall chemical composition is surprisingly stratified 
without significant mixing among Fe, Si and other ambient ejecta gas. Light echo (Rest 
et al. 2005) and X-ray observations (Hughes et al. 1995) yielded estimates of the age of 
N103B of ~ 850 and ~ 1000 - 2000 yrs, respectively. Given the age, the clump-remnant 
interaction should have started in N103B. The consistency between the ionization time scale 
ratio of Si to Fe and the density contrast across the Ni bubble shell demonstrated that 
the Ni bubble expansion had occurred, and yet the globally stratified composition reveals 
suppression of mixing in the clump-remnant interaction. This outcome is further consistent 
with our expectations of the exponential model. 

6. SUMMARY AND CONCLUSIONS 

We have investigated the structure and expansion properties of the Ni bubble shell 
in SNe la due to radioactive heating caused by ^^Ni ^^Co ^ ^^Fe decay sequence. 
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assuming an exponentially-declining density profile of the ejecta substrate. Based on the 
indication that ejecta clumps arise from the breakup of components of the Ni bubble shell, 
whether the inferred shell properties are compatible with those revealed for the ejecta knots 
in Tycho's remnant was studied. Adiabatic solutions were first obtained, and then radiative 
solutions, including the radiation diffusion process across the Ni bubble shell, up to the 
approximate frozen-in age ~ 10® s, were applied to estimate the clump properties. Since a 
realistic elemental distribution of ^^Ni results in a higher expansion velocity than a centralized 
distribution of ^^Ni used in our model, which effect is equivalent to the use of a lower initial 
^^Ni density contrast in the range of 0.5 < a; < 1, our u; — 1 cases may only hint at the 
lower limit on the clump speeds. We believe that our inferred clump strength is at least not 
overestimated. 

The expansion of the Ni bubble sweeps up a dense thin shell of ~ 0.1 Mq shocked gas in 
a typical Type la SNR. The density of the shell increases inwardly toward the bubble-shell 
interface, such that the highest computed density is limited by numerical resolution. Since 
the exponential profile is close to power-law profile that evolves from a steep power law to 
a fiatter one, the adiabatic solutions of the shell were approximated as self-similar solutions 
for a pulsar bubble shell in power-law ejecta of power index n with the shock propagated at 
an expansion rate a, which in our intermediate-class models falls in the range n = 2 — 3 and 
a < 1.07. The shell is thicker than given by the flat ejecta model, representing < 1% of the 
forward shock radius and the mean density contrast is enhanced by 50%. In the radiative 
case, the shell is broader and less dense; the total swept-up mass of the shell is ~ 50% less 
without the radiative precursor, but the shell's expansion velocity is only < 3% lower close 
to the frozen-in age. The radiative diffusion process does not substantially affect the clump 
speed or the initiative age of the clump-remnant interaction. 

A wide range of explosion parameters that are similar to those used in valid 1-D explosion 
models was explored. The exponential model gives a higher shell expansion rate and thickness 
ratio than the uniform model, and a lower shell velocity and swept-up mass. The expansion 
of the shell is only slightly faster than the free expansion of the ambient ejecta substrate. 
Consequently, the frozen-in velocity of the shell often rises with the velocity scale height Ve 
of the exponential profile - which characterizes the free motion of the ejecta substrate - or 
drops with the ambient ejecta density gradient. The velocity and swept-up mass of the shell 
appear to be inversely correlated: steeper models (with smaller velocity scale heights) tend 
to acquire a lower expansion velocity and a higher swept-up mass. 

Under intermediate explosion conditions with an ambient ISM density in the range of 
nam = 0.5 — 1 cm~^, the clump-remnant impact is anticipated to occur at ~ 300 yrs after 
the supernova is formed. Although the impact precedes Tycho's present state, the clump 
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speed (~ 6000 km for u; — 1) may not be sufficient. A higher clump speed and an earher 
clump-remnant impact is expected when the model incorporates a more realistic elemental 
distribution. The explosion parameters obtained from the successful explosion models such 
as the deflagration W7 and the delay detonation DD200 yield the most favorable result for 
Tycho's knots. A sub-Chandrasekhar-mass, low-cncrgy scenario such as the He detonation 
HeD6 may also satisfy Tycho's condition. A Chandrasekhar-mass and low-energy model 
such as the PDD is unfeasible, indicating that if Tycho's SN was intrinsically under luminous, 
then a near-Chandrasekhar-mass progenitor could not arise. Likewise, SN 1006 could not 
originate in a W7-like scenario since the currently-expected clump-forward shock impact is 
not observed. Recent studies nonetheless identified both Tycho's SN and SN 1006 as normal 
Type la supernovae. 

The global approximations of Type la models with exponential models may not be 
sufficiently accurate enough close to the contact discontinuity, particularly for the PDDlc 
model, where the ejecta appears to diverge from the exponential model at low velocities. If 
the ejecta density profile deviates from the exponential or power-law models, a higher clump 
speed and an earlier clump-remnant impact is expected for a flatter ejecta substrate, and vice 
versa. However, signiflcant deviations from the exponential model are not anticipated for the 
W7 and DD200c models. The clump-remnant impact is estimated to occur within ~ 2000 
yrs under all plausible explosion conditions of SNe la. As clumps are expected to be present 
only in a thin shell at the boundary of the Ni bubble, not throughout the bulk of the ejecta, 
the small clump size and a flat bubble density contrast should not induce vigorous mixing in 
the ensuing clump/bubble-remnant interactions, which such characteristic distinguishes the 
remnants of SNe la from those of the core-collapse SNe. 
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Fig. 1. — Evolution of the characteristic velocities of an SNR that arises from interaction of 
exponentially-declining ejecta with a uniform ambient medium, as described in DC98 and 
WCOl. Dotted line: Flow velocity of ejecta immediately below reverse shock (RS). Solid 
lines: Pattern velocities and deceleration parameters S — dlnr/dlnt (inset) of forward shock 
(FS), contact discontinuity (CD) and reverse shock. 
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Fig. 2. — Distributions of material density, velocity and thermal or radiative energy density 
of Ni bubble in standard exponential model at 5 x 10^ and 10^ s. Solid and dash-dotted lines 
represent adiabatic and radiative diffusive cases, respectively. Also plotted as dotted line is 
unaltered exponential density distribution. HD runs exploit a grid of 16000 uniform zones, 
resolving the shell into > 100 zones, while RHD runs exploit a grid of 8000 uniform zones. 
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Fig. 3. — Evolution of properties of Ni bubble shell in various models. The grids use 
16000 uniform zones for all runs, (a) Radius Rsh evaluated at shock front, (b) Expansion 
rate a = dlnRgh/dlnt. (c) Thickness ratio j3. (d) Density contrast (middle) of bubble 
interior near contact discontinuity, relative to the ambient ejecta substrate. In this plot, 
the PDD3, W7 and DD200c models appear to yield matching results, (e) Flow velocity 
Vfiow,sh = Rsh/t evaluated at densest point of shell (CD), (f) Swept-up mass M^. (g) Surface 
density as = Ms/AttRI,^. (h) Power index n of the density of the ejecta substrate at contact 
discontinuity. 
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Fig. 4. — Ni bubble density profiles at 10^ s (beginning of simulations, top and right axes) 
and 10^ s (bottom and left axes) for the exponential u; — 0.7 W7 model (solid hne), RHD 
counterpart (dashed-dotted line), and modified exponential model using initial ^^Ni mass 
fraction given by Nomoto et al. (dotted line). The RHD model uses grid with 2000 uniform 
zones and HD models use grid with 4000 uniform zones. 
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Fig. 5. — Ni bubble density profiles in RHD models at their estimated frozen- in ages. Models 
from left to right are, respectively, HeDlO (2 x 10^ s), HeD6 (7 x 10^ s), W7-2 (1 x 10^ s), W7 
(1 X 10^ s) and DET2env6 (1 x 10^ s). Sohd lines represent adiabatic cases and dash-dotted 
lines represent radiative cases. Model HeDlO uses grid with 2000 uniform zones while all 
other runs use grid with 4000 uniform zones. 
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Fig. 6. — Top: Flow velocities Vfiow,sh = Rsh/t evaluated at CD in radiative models. Models 
from top to bottom are the HeDlO, u; = 0.7 W7 (both HD and RHD cases), W7, DET2env6, 
HeD6, std (HD and RHD cases) and W7-2, respectively. 
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Fig. 7. — Hydrodynamical solutions of Ni bubble shell in standard model at 2.0 x 10^ s over- 
plotted on self-similar solution for n — 2.448 and I — 0.8958 (corresponding to a = 1.041). 
These self-similar solutions were obtained by integrating the self-similar differential equations 
using the fourth-order Runge-Kutta method. The shell is resolved into ~ 140 grid zones in 
the HD solutions and 200 zones in the self-similar solutions. The HD solutions exhibit large 
oscillations in the velocity and pressure distributions behind the CD, which are caused by 
the use of the shock-capturing scheme which aims to resolve an abrupt density change in 
a molecular length scale into a few grid zones, whose numerical origin is the same as that 
of the postshock oscillations. Traditionally, weaker shocks and higher numerical resolution 
are thought to produce more severe postshock oscillations, whose presence nevertheless de- 
pends on the problem; satisfactory design criteria that ensure that the captured shocks are 
simultaneously narrow and free from oscillations do not exist. For a discussion of postshock 
oscillations, see Arora and Roe (1997). 
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Fig. 8. — Variations in self-similar thickness ratio and mean density contrast with expansion 
rates for four density power indices. The lines from top to bottom refer to cases of n = 3, 
n — 2.4, n — 2 and n — 0. The solutions were obtained using a fixed integration step size of 
6.6 X 10~^. (Radius of shock front is one). 



-35- 





0.04 








8 


N 


0.03 






1 : 


6 














s 


0.02 






/' ; 
















G 














0.01 








2 









- .- -1 - .- : -.-.-\-.-.-. : r r. : .t r. : .j .r.n .-.t\.-. 








0.1 0.5 1 2 3 

^'rs 

Fig. 9. — Initial clump size (solid and dotted lines) and forward shock impact time (triple- 
dot-dashed line) as functions of reverse shock impact time. Solid lines: Size as fraction of 
intershock width. Dotted lines: Size as fraction of forward shock radius. Lines from top to 
bottom refer to cases /5 = 0.007, 0.006, 0.005, 0.004, 0.003, 0.002 and 0.001, respectively. 
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Fig. 10. — Cloud crushing time (solid hne) and ratio of cloud-crushing time to intershock- 
crossing time (dotted line) as functions of reverse shock impact time, assuming tec — 
f^fRsX^^'^/i'^f - vrs), P = 0.04 and x = 100 for the clump. 
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Table 1: SN la EXPONENTIAL MODELS 
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Quantities given in columns are as follows: (1) name of model; names with * are cases with ui = 0.5 (2) 

supernova mass (Mq) (3) supernova kinetic energy (lO^'^ergs) (4) initial mass of ^^Ni (Mq) (5) constant of 
the exponential model (10^ cgs units) (6) velocity scale height (km s^^) (7) frozen-in velocity (km s~^) (8) 
normalized frozen-in velocity (9) normalized age of the reverse shock impact (10) age (yr) of the reverse shock 
impact for an ambient density of Uam = 1 and 0.5, respectively (11) normalized age of Tycho's remnant for 
an ambient density of nam = 1 and 0.5, respectively. 



